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Introduction 

There are many systems which are described by nonlinear operators. A quadratic stochastic operator 
(QSO) is one of the simplest nonlinear cases. A QSO has meaning of a population evolution operator and it 
was first introduced by Bernstein in [1]. For more than 80 years, the theory of QSOs has been developed and 
many papers were published (see e.g. [4]-[7]). In recent years it has again become of interest in connection 
with its numerous applications in many branches of mathematics, biology and physics. 


Let E= {1, Deca m} be a finite set and the set of all probability distribution on E 


i=l 


SS ix = (2,,%,,.-.%,,) €0 "x, 2 0, for any i and Dx, = ' 
(1) 
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_ . Om m-1 
be the (m eer simplex. A QSO is a mapping defined as VS" 3S of the simplex 


eee m-1 
into itself, of the form V(x) =x eS , where 
i=) Foxit LEE, 

7 (2) 

= te 

and the coefficients * satisfy 
= >0O =] 
ijk jik — ~? ijk o 
ae forall &J © E. (3) 


(n) 
x (0) m-1 
The trajectory (orbit) - , of V for an initial value X © S is defined by 
SV eave). wa Ol as 


One of the main problems in mathematical biology is to study the asymptotic behavior of the trajectories. 
This problem was solved completely for the Volterra QSO. 


keé{i,j\, i j,.keEk 


The operator V is called Volterra QSO, if Fig =) for any . For the Volterra 


QSO the general formula was given in [4], 


' — 
Xx, =X, [ i Yous 
i=l 


; (4) 
a, =2P 


where “ ikke 


1 |< 


| for all kek 


ey and 4 


for b#K and Ge = 0. Moreover, “ki 

In [4], the theory of Volterra QSO was developed using theory of the Lyapunov functions and 
tournaments. But non-Volterra QSOs were not completely studied. Because, there is no general theory that 
can be applied for study of non-Volterra operators. 

In [7] Separable Quadratic Stochastic Operators (SQSOs) where introduced. Volterra QSO (4) has a 
form as SQSO, but in [8] it was proved that it coincides with a SQSO if and only if it is a linear operator. 

In recent years, Cubic Stochastic Operators (CSOs) have begun to be studied, which different from 
quadratic operators. 

In this paper we consider another class of cubic operators which we call separable cubic stochastic 
operators. 

In Section 2, we recall the definition of CSOs and definitions and known results. In Section 3 for a SCSO 
defined on the two-dimensional simplex, we prove that it has three fixed points and we find conditions on 


2 
parameter under which a fixed point is a repelling, attracting, or saddle point, the boundary OS” is an invariant 
set. 
Preliminaries and known results 
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Separable quadratic stochastic operator. Let us recall some necessary definition and notations. In [7] 
Separable Quadratic Stochastic Operators were introduced as follows: The QSO (2), (3) with additional 
condition 


=a,b ijJkKEE 


Fi ik jk for all 


(5) 


0, €U =(a, B=(b. - 
where Ay Dy =o entries of matrices A (4, ) and ( ik such that the conditions (3) are 
satisfied for the coefficients (5). 


Then the QSO V corresponding to the coefficients (5) has the form. 


x, =(V(x)), =(A(x)), (B (x), (6) 


where i=l ; jt 
Definition 1. [7] The QSO (6) is called separable quadratic stochastic operator (SQSO). 


m 


>i, El 


oa es as roy. ae . 
From the conditions '* and <=! for all a, it follows that the condition on matrices A 


a,b. >0 T T 
and B that * J ; AB = 1, where B is the transpose of B and 1 is the matrix with all entries 
OZ . Bab, b. ) 
1’s. If ° (2...) is the ! -th row of the matrix A and iim) is the J -th row of 
T 
the matrix B , then from AB’ =1 ve get 


ap =1 for all i, j=l,...,m 


Hs 


For a fixed J , the above condition implies that 


A(b!) =(1,1,...1) 


(7) 
. 
If det( A) - 0 then (7) gives pe = b| for all kJ E ie., all rows of B are the same, 
therefore det(B) =U Similarly, if det(B) ia 0 then all the rows of A must be the same, so 
det(A)=0 


‘. m-| m-| 
Cubic stochastic operator. The CSO is a mapping Ws > Sof the form 
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i,j,k=1 (8) 
where ‘** are coefficients of heredity such that 


SY et eS Le 8 es ee ee BO 
i (9) 
i, jk 


P. 
and we suppose that the coefficients /* do not change for any permutation of 


Note that W is a non-linear operator, and its dimension increases with mM. Higher-dimensional 
dynamical systems are essential, but only relatively few dynamical systems have yet been analysed. 


(n) 
ee Oe Gm RE ty | | 
For a given , the trajectory n20 of initial point X ~ under action of CSO (8) is 


x -w ey ) (0) o(x” ) 
, where n=0,1,2,... with X =X. Denote by the set of limit 


{x” \ x” \ C St dy 
points of the trajectory n=0 , Since n=0 and S is a compact set, it is follows that 


o(x'”) #D = o(x') 


defined by 


foe) 


(0) ) 
oo(x 
consists of a single point, they the trajectory converges and is a fixed 


m-I x)=x Fix 
point of the operator Wia point X¥ © Ss” is called a fixed of the W ig w( ) . Denote by (W) 
the set of all fixed points of the operator W , Le. 


Fix(W) = {x eS"! :W(x)= x} 


DW (x’) = (OW, / Ox, )(x') 


Let be a Jacobian of W at the point x 


ae 2 — _ DW(x’) | 
Definition 2 ([3]): A fixed point X is called hyperbolic if its Jacobian has no eigenvalues 


on the unit circle in 


Definition 3 ([3]): A hyperbolic fixed point x’ is called: 
Pe |  DW(x’) 
(i) attracting if all the eigenvalues of the Jacobian are in the unit disk; 


D (x" 
(ii) repelling if all the eigenvalues of the Jacobian are outside the closed unit disk; 


(iii) a saddle otherwise; 


© 2023, CAJMTCS | CENTRAL ASIAN STUDIES www.centralasianstudies.org ISSN: 2660-5309 | 84 


CENTRAL ASIAN JOURNAL OF MATHEMATICAL THEORY AND COMPUTER SCIENCES _ Vol: 04 Issue: 12 | Nov 2023 


Main result 


In this section we consider CSO (8), (9) with additional condition 


ae = a,b jl Cy 
? 


i,j KLEE 


for all (11) 


lu © A=(a,) B=(b,) ad C= 


a,,b . 
where ”’ entries of matrices 
conditions (9) are satisfied for the coefficients (11). 


(Cy ) such that the 


Then the CSO W corresponding to the coefficients (11) has the form 


x! =((x)), =(A(x)), (B(x), (C(%)) o 
(A(x)), = as (B(x)), _ m DX: (C(x)), = Meu 


where ae and : 


Definition 5. The CSO (12) is called separable cubic stochastic operator (SCSO). 


Lemma 1. The condition (11) is sufficient for a CSO to be product of three linear operators, but the 
condition is not necessary. 


P A=(a,) B=(b,) 
The CSO with coefficients “*-’ can be written as the product of three matrices (4, il 


l,j KLEE 


and C= (Cy ) if and only if, for any it holds 


desea se tT Ae Di PP Py = G0 Cy G0, C+ G,0,C , + 


ijkl kij. ikj, Kji.l jik,l jkil il’? jl kl Kil i? kL jl 
GON, oa a b.Cy as GD, (13) 


Thus the condition (11) is a particular case of (13). 


Let us consider the following matrices: 


100 a, 2 2S i tee 4 
Aa=!6 f OL. Bealia 1 t es) 7 1 1+al, 
0 0 1 1 l-a 1 1 a 1 1 
2 (18) 
where ae[-L]] 


.@2 2 
Then corresponding SCSO W:S° >S* is. 


© 2023, CAJMTCS | CENTRAL ASIAN STUDIES www.centralasianstudies.org ISSN: 2660-5309 | 85 


CENTRAL ASIAN JOURNAL OF MATHEMATICAL THEORY AND COMPUTER SCIENCES _ Vol: 04 Issue: 12 | Nov 2023 


Ba =e +(1—-a)x, +x)[ +X, +(1-$}s} 


W 34x) = +X, +(1-a)x,)((l+a)x, +X, +x), 


xi = «(1+ ae +X, + ](s +(1+a)x, +x,). 
2 

(19) 

X, +X, +x, =1 


Using the equation we rewrite the operator (19) as follows 


(20) 


Evidently, that if ¢ = 0 the SCSO (20) is the identity map. For this in the below, we consider the case 
when @ #0. 


Tl, ={xeS?:x,=0,i¢ac{l,2,3}} 


a 


2 
Let a face of the simplex S be the set 


int S* ={xeS’: xxx, >0 2 eas 
as oan let the set os” = S* \int S be the interior and 


e, =(1,0,0) e, =(0,1,0) e, =(0,0,1) 


Let the set 


2 
the boundary of the simplex S , respectively. Let 
vertexes of the two-dimensional simplex. 


be the 


Theorem 1. For the SCSO W (20), the following assertions true: 


| eran By 


2 
(i) The face {1,2} ‘ {1.3} ; ee) of the simplex S are invariant sets; 


Fix(W) = {e,,e,,e,} 


(11) 
un pA €|—-L0 ; ; . ; ; ; 
(iii) If | : ) , then ©) is an attracting point, © is a saddle point and © is a repelling point; 
aeé(0,1 
If ( ° , then <i is a repelling point, © is a saddle point and ae is an attracting point. 


Proof: (1) Obviously. 
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W(x)=x 
(ii) To find the fixed points we consider the equation ( ) , that is the following system of 
equations 


(21) 
(a) Let “= 0. Then from the second equation of (21) it follows 
X, =X, —AX,X, = x, =0 and 1=1—-ax, 
It is clear that if a 0 then we have = | If a >0 


then from the last equation one has 
x,=0 = x, 1 


We take solution a= 1 then it follows that = 0 Consequently, if 4 0 we obtain the fixed 
points ae ce 


? 


Similarly, in the case xy = 0 and ds 0 we have the fixed points it & and in ) respectively. 


(b) Suppose that Piyhahe 0. Then from the system (21), one has 


(22) 


aé|-1,0 : 
If | ) then from the first equation of (22) we get 


a 1 
I=(1-ax)(1-$ he aH +—x 


=> 2 3 
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E —] = = 
But from " | .0) and from the last equation it follows that ae 0 and a, 0 this contradicts 
to 21% FO. 
€(0,1 
a (0, , then from the third equation of (22) we have 
a a 1 
1=| 1+—x (1+ ax, ) —XxX,=-| —x,+x 
9 1 2 9 17%2 9 1 2 
=> 
E — = 
But from e (0, and from the last equation it follows that ae 0 and My 0 this contradicts to 
X,X,X, #0 


Fix(W) aint S? =@ 


Consequently, we have that 


(ii1) To find the type of fixed point of the SCSO (20), we rewrite it in the form 


Se (1max,)(1-S-% -x)}, 


x, = x,(1+ax,)(1-a(1-x, —x,)). (23) 


(x,,x,)e{(x,y):x%,y20,0<x+yslh 


X,,X. : . 
where and ~1!’° 2 are the first two coordinates of a 


: ee : S 2 
point lying in the simplex : 


a 
P A, =1+=— 
The Jacobian of the operator (23) at a fixed point ~ has the following eigenvalues 2 ; 


A, =l+a at the ~ has the following eigenvalues A, - l-a A, =I+a and at the . has the 
a 
Se 
following eigenvalues “1 ; 2. Therefore, it follows that: 
€|—1,0 
If ° : ) then, the fixed point “i is a attracting point, ve is a saddle fixed point and the fixed 


point Es is a repelling point. 


If = 0 then, fixed points “i. © and a are non-hyperbolic points. 
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If 


a <(0,1| 


then, the fixed point “i is a repelling point, “ is a saddle fixed point and the fixed point 


c.. : ; 
3 is a attracting point. 


The theorem is proved. 
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